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The relative roles of divergence and velocity slip in the stability 
of plane channel flow 
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PACS 47 . 20 . -k - Flow instabilities 

PACS 47. 15.Fe - Stability of laminar flows 

PACS 47 . 27 . nd - Channel flow 

Abstract. - Wall slip and wall divergence are known to have large and opposing effects on the sta- 
bility of flow in a two-dimensional channel. While divergence hugely destabilises, slip dramatically 
stabilizes the linear mode. In a non-parallel stability analysis, we study a combination of these 
two effects, since both will coexist in small-scale flows with wall roughness. Our main results are 
(i) that the stabilising effect of slip is reversed at higher angles of divergence, (ii) transient growth 
of disturbances is unaffected by either wall-divergence, or by slip at any divergence. Moreover, 
at the Reynolds numbers relevant here, transient growth is too low to be a significant player in 
transition to turbulence, which is more likely to be driven by linear instability. 



Introduction. — The laminar flow through a two- 
dimensional channel is linearly stable up to a Reynolds 
number of 5772, based on the channel half- width and the 
Centreline velocity. Since the flow usually goes to turbu- 
lence at a much lower Reynolds number (between 1000 
and 2000), the mechanisms leading to turbulence are ac- 
cepted to be either directly nonlinear, or triggered initially 
by transient algebraic growth due to the superposition of 
linearly decaying waves. In the latter case, nonlinearities 
are said to take over once transient growth increases dis- 
turbance amplitudes to significant levels. 

In this paper we investigate the combination of two large 
and opposing contributors to stability; channel divergence 
and velocity slip. We study both exponential and transient 
growth of instabilities. Our motivation is that in small- 
scale flows, wall-roughness can create local divergence and 
convergence of the flow, and slip can be signiflcant, espe- 
cially if the walls are hydrophobic, so the two effects are 
likely to act simultaneously. It was noted decades ago [1,2] 
that a small divergence of the wall can cause the critical 
Reynolds number for linear instability to plummet by two 
orders of magnitude. The transient growth of disturbances 
in this flow has not been investigated before, to our knowl- 
edge. The instabiliy of channel flow in the presence of ve- 
locity slip at the wall is receiving attention recently ( [3,4] 
and references therein) . The consensus for a plane channel 



is that wall slip hugely stabilizes the linear mode, but does 
not affect transient growth of instabilities much [5-7]. 

Under no-slip at the wall, our non-parallel linear stabil- 
ity analysis conflrms the result obtained by [1] under the 
parallel flow assumption, that with increase in divergence, 
the critical Reynolds numbers falls sharply. There is a 
small quantitative change in the result due to nonparallel 
effects. On the other hand, we show that for small lev- 
els of wall divergence, transient growth of disturbances is 
unaffected by wall divergence. At the relevant Reynolds 
numbers, the magnitude of maximum of transient growth 
is small. Slip at the walls does not in this case have the 
"obvious" effect of huge stabilisation. In fact, wall slip can 
further destabilise the flow, although slightly. 

At low Knudsen numbers, we take the Navier-Stokes 
equations to be valid in conjunction with the slip boundary 
condition of Maxwell [8] , 



U±Kn^ 

ay 



0, a.t y = ±1. 



(1) 



This model is correct to 0(1) in the Knudsen number, Kn, 
expressed here as the ratio of slip length (Ig) to the local 
half- width II{x) of the channel. The coordinates along, 
and normal to, the channel centreline are denoted as x and 
y respectively. All quantities are non-dimensionalised us- 
ing H{x) and the centreline velocity Uc{x). In particular. 
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Fig. 1: (a) A schematic diagram showing the coordinate sys- 
tem, slip length and slip velocity, (b) Variation of U" at the 
wall for different S (solid line: Kn = 0.1 and dashed line: 
Kn = 0). 
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Fig. 2: (a) Streamwise, and (b) normal, components of the ve- 
locity at various values of the divergence parameter 5*. Solid 
line: Kn = 0.1, dashed line: no slip. The U and V are sym- 
metric and anti-symmetric about the channel centreline {y = 0) 
respectively. 



the relationship da;^ = Hdx is used later, where the sub- 
script d represents a dimensional quantity. The Maxwell 
slip model breaks down around Kn ~ 0.15 [9], so we, like 
others [5-7, 10] who use this model, restrict ourselves to 
values below this. The slip length and slip velocity {Us) 
are shown schematically in figure 1(a). 

Incidentally the work of Lauga & Cossu [5] set to rest 
a recent controversy over the effect of slip on the stability 
of channel flow. Chu [10-12] had found a large destabil- 
isation due to slip, in contrast to the earlier prediction 
of Gersting [13] of a huge stabilisation. The conclusions 
of [13] are confirmed by [5], who showed that there was 
a lapse in the application of the slip boundary condition 
in the work of Chu. There remains a small discrepancy 
in the effect of slip on transient growth, between the re- 
sults of [5] and [6]. Both agree that the effect of slip is 
small, but the direction of the effect is stabilising in [6] 
and destabilising in the other. We have checked that this 
discrepancy arises from a difference in the definition of 
Reynolds number, since the velocity scale (from their re- 
spective analytical expressions) used by [6] decreases with 
increasing slip while that of [5] goes up with slip. In the 
present paper we avoid this ambiguity by scaling with the 
centreline velocity, obtained numerically, as explained in 
the next section. For the sake of validating our results, 
when we employ the respective definitions of [5] and [6], 
our calculations for a plane channel agree very well with 
each. 

For a divergent channel, we use the Jeffery-Hamel mean 
flow profiles, described in the next section, followed in sub- 
sequent sections by the linear stability and computations 
of transient disturbance growth. 

The basic flow. — At low angles of divergence, the 
incompressible laminar flow in a wedge-shaped channel is 
given by the Jeffery-Hamel equation (see e.g. [1]): 

U'" + 2SUU' = 0, (2) 

which is valid for flows with slip boundary conditions as 
well. All effects up to the first order in the divergence 
have been retained in the Navier-Stokes equations, and 
similarity has been assumed. The primes denote derivative 
with respect to y. The only parameter is the product 



S = aRe, of the slope a of the wall and the Reynolds 
number Re. (Some of the results below are presented in 
terms of the angle of divergence 6 = tan~^ a). From global 
continuity we see that Re = H{x)Uc{x)/v is a constant for 
a given flow, v being the kinematic viscosity. The normal 
velocity 

V = ayU (3) 

from continuity. Note that V, being normal to the centre- 
line, is not zero at the wall when there is a non-zero slip 
(Us). Equation (3) with the boundary conditions given by 
equation (1) at the wall, and J7 = 1, ?7' = at the cen- 
terline, is solved by a fourth order Runge-Kutta method. 
Thus for different Knudsen number the profile is computed 
numerically. The effect of slip at different S on the stream- 
wise and normal components of the velocity is shown in 
figures 2(a) and (b) respectively. 

Ever since Rayleigh [14] proved that a point of infiexion 
in the velocity profile is a necessary condition for inviscid 
instability, the second derivative U" has been known to be 
a sensitive indicator of linear instability, especially when it 
goes through a zero somewhere in the flow. In a divergent 
channel, U" at the centreline remains negative, and for 
small divergences, is close to the value in a plane channel 
of —2. The dependence of U" at the wall on S is shown in 
figure 1(b). A positive value is indicative of the existence 
of an inflexion point. In the absence of slip, the profile is 
inflexional beyond S = 2.95, whereas for Kn — 0.1, the 
profile is infiexional for S > 2.5. Going by this inviscid 
argument, a divergent channel might be expected to be 
more unstable with slip than without, but the difference 
is too small to predict the direction of the effect at finite 
Reynolds numbers. We may however estimate that for 
a profile that is near-infiexional, the divergence could be 
more important for stability than changes in the bound- 
ary conditions. In the next section, we will investigate if 
the above argument is correct with a non-parallel linear 
stability analysis. 

Non-parallel stability analysis. We consider 
small divergence angles and Reynolds numbers are much 
larger than 1 , so streamwise variations in the basic flow are 
slow. A consistent non-parallel approximation is to retain 
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all terms up to 0{a) and O (i?e^^) and to neglect terms 
of higher order in either of these. As in parallel stability 
analysis [15], each flow quantity is expressed as the sum of 
a steady mean and a time-dependent perturbation, such 
as 

u=^U{y) +u{x,y,z,t). (4) 

Note that, y being the similarity variable, the basic ve- 
locity in equation 4 is a function of y alone. Our nu- 
merical results for non-similar velocity profiles show that 
U ~ U{y) is a. fair assumption whenever dU/dx is smaller 
than the inverse of the Reynolds number. The depen- 
dence on the streamwise coordinate may be expressed as 
the product of a rapidly varying wave-like part scaled by 
a relatively slowly varying function (see e.g. [16, 17]), and 
normal modes may be used in the spanwise coordinate z 
and in time, as 



Real < u{x, y) exp 



a{x)dx + f3z — Lot 

" '(5) 

such that du/dx ^ O {Re The streamwise and span- 
wise wavenumbcrs are a and /? respectively and w is the 
disturbance frequency. The equations for linear stability 
under these assumptions are 

du . du dv , . 

-au + - — I- lau - ay— — \- - — h \pw = 0, (6) 
ox oy oy 



U 



du 
dx 



-\- {ia — 2a) u 



+ {v-ayu)U' (7) 



. . dp dp 
^2a~,a)p^-+ay- 



1 

Re 



dy^ 



dp 
dy 



U 
1 

'Re 



dv 
dx 



ayU'v 



(8) 



-iujw + U 



-i(3p' 



1 

Re 



dw 
dx 
^d^w 
dy'^ 



(9) 



The boundary conditions are 



dz 



dw 



u ± Kn— = V ^ au — w ± Kn—— — at y = ±1, 
oy dy 

(10) 

with finite pressure at the walls. In equations (6) to (10), if 
we set a, and all the streamwise derivative terms to zero, 
we obtain the Orr-Sommerfeld and Squire equations for 
parallel flow [15]. 
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Fig. 3: (a) Variation with velocity slip at the wall of the critical 
Reynolds number Rccr for linear instability in a plane chan- 
nel. Line: present study, symbols: [5]. (b) Neutral stability 
curve for different values of Kn, angle of divergence = 0.1°. 
(c) Neutral stability curve for different values of Kn, angle of 
divergence = 1°. (d) Re/Reo Vs Kn for different angle of 
divergence {6). 



We have followed two completely different approaches, 
described in detail in [17] and [18] respectively, to solve 
the stability problem of equations (6) to (10). The results 
agree very well for divergent flows without slip, but the 
latter method is not applicable for slip flows. In the flrst 
method the equations at two neighbouring streamwise lo- 
cations are rewritten in the form of an eigenvalue problem 
of larger dimension. This is made possible upon descretis- 
ing the streamwise derivative, which couples neighbouring 
streamwise locations to one another, and using the fact 
that the dimensional frequency remains constant down- 
stream. In the second approach, the lowest-order stability 
problem correct to 0{Re^^^^), comprising ordinary differ- 
ential equations in y is derived. This is solved first as an 
eigenvalue problem, exactly as one would solve the Orr- 
Sommerfeld equation. Higher order effects are obtained up 
to the desired order of accuracy by exploiting the proper- 
ties of adjoint differential operators. The disturbance am- 
plitude is expressible in terms of an ordinary differential 
equation in x, whose coefficients contain the lowest-order 
eigenfunctions. Both methods use Chebyshev collocation 
spectral discretisation. We always find that the even mode 
in V is much more unstable than the odd mode. When ap- 
propriate, this fact is used to speed up the computations 
by considering only a half-channel. 

While Squire's theorem is not applicable directly for 
non-parallel flows, it is usually the case that two- 
dimensional perturbations are the least stable in flows 
which vary slowly in x. This is true in the case of di- 
vergent channels too [19], so only two-dimensional pertur- 
bations are presented in this section on linear instability. 
The two parameters in the problem under consideration 
are channel divergence and wall slip. For a plane channel 
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Fig. 4: (a) Variation of the critical Reynolds number, Rccr with 
wall divergence, solid line: Kn = 0.1, dashed line: Kn = 0. (b) 
A portion of figure 4 (a) is magnified, with the corresponding 
results from the Orr-Sommerfeld equation shown in symbols. 



{6 = 0), the present solutions match well with [5] as seen 
in figure 3(a). It is seen that the critical Reynolds number 
increases sharply with increase in the Knudsen number 
(nearly 175% for a Knudsen number rise of 0.05). For dif- 
ferent wall divergence the same variation is given in next 
two figures as neutral stability curves. The boundaries of 
neutral stability at various Knudsen number for angles of 
divergence 9 = 0.1° and 1° are shown in figures 3(b) and 
3(c) respectively. At the lower level of divergence, slip is 
seen to stabilise the flow, but much less than it would in a 
plane channel, only ~ 22% for an angle of 0.1°, figure 3(b). 
At the larger divergence, increasing the Knudsen number 
has a small destabilising effect only ~ 10% for an angle of 
1°, figure 3(c). The effect of slip and divergence is sum- 
marised in the plot of the critical Reynolds number Rccr 
in figure 3(d). To highhght the result, we plot Recr/Rso 
Vs Kn, where Rcq is the critical Reynolds number for 
the no-slip {Kn = 0) case. As the angle of divergence 
increases, the effect of slip is progressively reversed, from 
being hugely stabilising to mildly destabilising. 

Figure 4 shows the cross-over in the critical Reynolds 
number. At a Knudsen number of 0.1, Rccr varies with 
wall slope approximately as a power law, with an exponent 
of —0.8. This means that the parameter 5* at which the 
fiow is neutrally stable increases slowly with the wall slope, 
roughly as S' ~ a"'^. It is worth noting that as S increases, 
the slope of the velocity profile at the wall decreases, so 
a given Knudsen number corresponds to smaller slip ve- 
locities at higher divergences. The results from the Orr- 
Sommerfeld equation applied to the Jeffery-Hamel profiles 
are shown in the same figure. It may be concluded that 
for wall divergences below a degree, fiow non-parallelism 
acts only through the basic profile and the explicit non- 
parallel terms are negligible. This has been checked to be 
true for the entire spectrum of eigenvalues for a variety of 
conditions. Beyond this angle, explicit non-parallel effects 
are noticeable but not large. Incidentally, this finding is 
in contrast to what is seen in pipe flows, where the non- 
parallel terms have a large destabilising effect [17]. At 
the largest divergences shown, i.e., for a ~ 0.2, the flow 
is mildly separated at the critical Reynolds number. Any 
flrm statements at these angles of divergence would re- 
quire a non-parallel stability analysis correct to a higher 




Fig. 5: The contour of Gmax for Re = 200 without slip (solid 
lines) and at Kn = 0.1 (dashed lines) for (a) a = 0.005, and 
(b) a = 0.02. At the latter wall divergence, a Reynolds number 
of 200 is 95% of the critical Re for linear instability. 



order of accuracy in a and O (i?e~^) than considered here. 
We only show Orr-Sommerfeld results as a qualitative in- 
dicator that nothing special happens when one encounters 
a mildly separated profile. 

Transient growth. — In contrast to the least stable 
linear cigcnmodc, the maximum transient growth of dis- 
turbance kinetic energy is shown here to depend primarily 
on the Reynolds number of the fiow and not on the slope 
of the walls. The introduction of velocity slip at the walls 
has a very small effect in the presence of divergence as well, 
similar to the findings for a plane channel of [5] and [6]. 
Given the negligible effect of the non-parallel terms on 
the stability, we use eigenvalues and eigenfunctions from 
the Orr-Sommerfeld and Squire equations with slip bound- 
ary conditions. If A = diag{L0i,u!2, ■■■Wn} is constructed 
out of the N least stable eigenvalues of this system, and 
K is similarly constructed from the expansion coefficients 
of the corresponding eigenfunctions, the resulting distur- 
bance kinetic energy, g{t), may be written as [20] 



k{0)\ 



(11) 



Maximising g{t) for all possible initial conditions k{0), we 
define 



G{t) = max g{t). 

K=iO 



(12) 



In a linearly stable situation, the quantity G{t) increases 
initially, and after attaining a maximum Gmax, decays to 
zero as i ^ CK). For a given value of Re, a, (3 and 6, Gmax 
is thus the maximum possible factor over which the ini- 
tial disturbance kinetic energy can grow algebraically. Our 
computations of G{t) for a plane channel at Reynolds num- 
bers of 1000 and 1500 (appropriately scaled) at different 
Knudsen numbers are, as mentioned above, in excellent 
agreement with [6] and [5] respectively. 

The contours of constant Gmax for a range of a and /3 
are plotted in figure 5 for a Reynolds number of 200 and 
wall slopes of 0.005 and 0.02. The solid fine represents 
the Gmax at Re = 200 for Kn = 0.0, whereas the dashed 
line represents Kn = 0.1. It is evident that velocity slip as 
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well as channel divergence have only a marginal effect. We 
now define a Reynolds number i?eo.95 equal to 95% of the 
linear instability-critical Reynolds number Rccr for that 
particular angle of divergence. The algebraic amplification 
of disturbances obtained at this Reynolds number is close 
to the maximum possible, since exponential growth will 
take over at a higher Reynolds number for that angle of 
divergence. The Reynolds number of 200 at which the 
contours are shown is about two-fifths of Rccr at a = 
0.005, while at a = 0.02, Reo.95 = 200. The point to be 
stressed here is that, even at this "high" Re for a = 0.02, 
the growth is comparable with that of a = 0.005. The 
maximum algebraic growth (i.e., the maximum contour 
value in figure 5) obtained just before linear instability 
serves only to amplify the initial perturbation by a factor 
of ~ 15. This is in contrast to a plane channel, where the 
factor is about 200 even at one-sixth of Rccr [20]. The 
results are qualitatively the same over a range of Re and 
0. 

Throughout the Knudscn number range considered 
here, the largest transient growth is obtained for stream- 
wise independent disturbances (a = 0) of spanwise 
wavenumber (3 close to 2, so the remaining results are 
presented at these wavenumbers. It should be stressed 
here that the comparisons made are valid irrespective of a 
and p. Figure 6(a) shows the Gmax a-t Reynolds numbers 
of 150 and 300. The slope at which the fiow becomes 
linearly unstable (to two-dimensional perturbations) is 
shown by the vertical line in each case. The values of 
Gmax{oi = 0,/? = 2) shown beyond this point correspond 
to the first local maximum. Although this value sharply 
increases with divergence, it has been checked to be ir- 
relevant in comparison to the exponential growth in the 
two-dimensional mode. In the linearly stable range, tran- 
sient growth is (i) small, and (ii) relatively insensitive to 
wall divergence, unlike the dominant linear mode, which 
is hugely destabilised. The small destabilising effect of 
slip on transient growth is summarised in figure 6(b). The 
combination of wall divergence and wall slip on plane chan- 
nel is, thus shown to have insignificant algebraic growth. 
We may conclude that the algebraic growth mechanism is 
not a major player in diverging channels, whether with or 
without slip. 

In figure 7(a), it is seen that the transient growth ob- 
tainable just below Rccr (i-e., /Jeo.gs) drops sharply with 
wall divergence, and is very small even for modest diver- 
gence. This result is compared to the maximum transient 
growth in a plane channel in figure 7(b), solid line. Note 
that the dashed lines in this figure consist of points from 
figure 7(a) obtained at various levels of wall divergence, 
each at its Reo.95- At values of wall divergence below 2 
degrees or so, the two lines without slip are almost indis- 
tinguishable, showing that Gmax is primarily dependent 
on the Reynolds number and not on the slope. Scaling 
arguments at small a show [20] that Gmax should vary 
as Re^. This has been shown to hold true for a vari- 
ety of fiows such as for plane Poiseuille, plane Couette, 
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Fig. 6: The maximum transient growth of disturbance kinetic 
energy at a = 0, /? = 2.0 as a function of (a) slope, and (b) 
Knudsen number. The Reynolds number in (b) is 200. The 
vertical lines in (a) show the divergence beyond which the flow 
is linearly unstable at each Reynolds number. 



flow in circular pipes and Blasius boundary layers, with 
different constants of proportionality. The relationship is 
seen to hold in the case of divergent flows as well, with 
Gmax ^ f{Kn)Re^ where /{Kn) is a weak function of 
Kn. At higher divergences, where the velocity profile is 
close to separation at the wall, there is a clear deviation 
from the power law. This could be of theoretical interest, 
and merits further investigation. 

For the streamwise-independent modes {a ~ 0) our find- 
ing that the velocity profile, and thus the level of diver- 
gence, do not affect the transient growth, are consistent 
with earlier analyses on boundary layers [20,21]. How- 
ever, at this point we do not have an analytical expla- 
nation for the insensitivity to divergence. Note that the 
Orr-Sommerfeld and Squire equations reduce in this case 
to 



1 



Re \ dy^ ^ dy"^ + P v 



t^- I TT^ - P V 



1 

Re 



P^v- 



dy^ 



ApU'v. 



(13) 



(14) 



fj is the y component of disturbance vorticity whose am- 
plitude 77 is defined as in equation (5), and uj = iwi, since 
ujr ^ for a = 0. v is now independent of the velocity 
profile, but rj is not. 

For Kn > 0.1 wc find that the nature, and level of 
stabilisation is different, but since the slip model is no 
longer appropriate here, this range is being approached in 
a different study. 

To summarise, it is known that in a plane channel, at 
Reynolds numbers of about a quarter of Recr, the levels of 
transient growth are so high that they can trigger nonlin- 
earities and a transition to turbulence [21]. Recent work 
shows that a velocity slip at the wall (for small Knudsen 
numbers) does not much affect transient growth. The im- 
plication is that, were the channel to be perfectly plane, 
one would not expect any big difference in the route to 
transition, with or without slip. Here, we have investi- 
gated how two opposing sources of linear instabilty affect 
channel flow. We show that for small divergences of the 
channel wall, the transient algebraic growth is an insignifi- 
cant player, whether or not the flow is slipping at the wall. 
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Fig. 7: Gmax at Of = and P — 2, (a) as a function of wall 
divergence at Jieo.95. The point at the lowest slope (a = 10"'') 
is obtained at a Reynolds number of -Reo.95 = 4282, while at 
a = 0.2, -Reo.95 = 56.6. (b) as a function of Reynolds num- 
ber. The solid line is for a plane channel without slip. The 
long dashes are the quantity shown in (a) plotted now against 
Reo.95- The coincidence of the two lines over most of the range 
shows that unlike linear instability, the transient growth is not 
affected by wall slope. The power law behaviour is seen to hold 
for slip flow as well. 

Linear instability on the other hand occurs at Reynolds 
number two orders of magnitude lower than in a plane 
channel in spite of the presence of wall slip. We predict 
that in the presence of divergence and slip, route to chaos 
is via linear instability for channel and hope to motivate 
experimental verification of this assertion. Further, we 
predict that in the presence of local divergence a finite 
amplitude growth of disturbance will lead to turbulence, 
irrespective of background noise or slip. 

We have considered only Jeffery-Hamel flows, represent- 
ing long channels with constant divergences, but our re- 
sults are relevant to channels of arbitrary geometry with 
local divergences extending to a few channel widths. We 
have checked that numerical profiles obtained in such 
channels match closely with Jeffery-Hamel profiles. Re- 
stricting ourselves to situations where the inlet flow is rel- 
atively quiet, we have not touched upon the role of non- 
linearity here, but there are several interesting questions, 
such as whether a nonlinear self-sustaining mechanism of 
the general class of [22] would be in operation or whether 
nonlinear traveling-wave solutions [23,24] will be observed. 

* * * 
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